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Abstract 

A  two-dimensional,  isothermal  mathematical  model  of  an  H2-C12  single  fuel  cell  with  an  aqueous  HC1  electrolyte  is  presented.  The  model 
focuses  on  the  electrode  reactions  in  the  chlorine  cathode  and  also  includes  the  mass  and  momentum  balances  for  the  electrolyte  and  cathode 
gas  diffusion  layer.  There  is  good  agreement  between  the  model  predictions  and  experimental  results.  Distributions  of  physical  parameters 
such  as  reactant  and  product  concentrations,  solution  and  solid  phase  potentials  and  local  current  densities  and  overpotentials  as  a  function 
of  cell  voltage  are  presented.  Effects  of  varying  the  initial  electrolyte  concentration  and  operating  pressure  are  analysed.  It  was  found  that  an 
electrolyte  inlet  concentration  of  6  mol  dm-3  gave  the  best  cell  performance  and  that  an  increase  of  operating  pressure  gave  a  steady  increase 
of  the  fuel  cell  performance. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Fuel  cells  are  in  general  used  as  generators  of  electrical 
energy,  either  for  mobile  or  stationary  applications.  In  a  H2- 
CI2  fuel  cell,  oxygen  is  replaced  by  chlorine  as  the  oxidiz¬ 
ing  agent;  hence,  the  product  from  this  cell  is  not  water,  but 
hydrogen  chloride.  Consequently,  the  H2-CI2  fuel  cell  is  as 
much  an  electrochemical  reactor  as  a  generator  of  electrical 
power.  The  use  of  H2-CI2  fuel  cells  for  co-production  of  hy¬ 
drochloric  acid  and  electrical  power  in  industrial  applications 
has  been  proposed  for  plants  having  an  excess  of  hydrogen 
and  chlorine  readily  available;  the  chlor- alkali  industry  has 
a  potential  for  energy  savings  with  the  use  of  H2-CI2  fuel 
cells  [1].  The  use  of  H2-CI2  fuel  cells  in  industrial  processes 
where  chlorine  is  produced  as  a  by-product,  for  example,  in 
magnesium  electrolysis,  can  also  be  beneficial.  H2-CI2  fuel 
cells  have  also  been  proposed  for  space  applications  [2]  and 
for  distributed  energy  storage  systems  [3-5]. 
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A  mathematical  model  of  a  fuel  cell  provides  a  better  un¬ 
derstanding  of  the  complex,  coupled  phenomena  that  occur 
in  such  systems,  as  well  as  useful  information  for  scale-up 
and  design.  It  also  enables  prediction  of  the  cell  performance 
as  a  function  of  operating  conditions.  To  our  knowledge,  no 
mathematical  models  of  H2-CI2  fuel  cells  have  previously 
been  presented  in  the  literature.  However,  the  system  has 
strong  resemblances  to  the  alkaline  fuel  cell,  of  which  sev¬ 
eral  models  have  been  presented  [6,7].  In  addition,  the  system 
also  has  similarities  with  the  PEM  fuel  cell  which  has  been 
extensively  modeled  [8-17]. 

In  this  paper,  a  two-dimensional,  steady-state ,  isothermal 
model  for  a  single  cell  of  a  hydrogen  chlorine  fuel  cell  is  pre¬ 
sented.  The  model  considers  four  of  the  six  layers  of  the  H2- 
CI2  single  cell:  membrane,  separator,  cathode  catalyst  layer 
and  the  cathode  gas  diffusion  layer.  The  mathematical  model 
is  solved  using  the  Femlab  3.0a®  program  package  and  the 
chemical  engineering  toolbox  from  Comsol  AB.  Profiles  of 
local  overpotential,  current  density  and  electrolyte  concentra¬ 
tion  are  obtained  as  a  function  of  cell  voltage  and  electrolyte 
concentration. 
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Nomenclature 

ag  specific  area  of  the  gas-electrolyte  interface 
(cm2  cm-3) 

af  activity  of  species  i,  adsorbed  on  electrode  sur¬ 
face 

ctf  activity  of  species  i,  in  bulk  solution 
Acat  catalyst  surface  area  (cm2  g-1) 

Cq  electrolyte  concentration  (mol  cm-3) 

cg  molar  concentration  of  an  ideal  gas  at  a  specific 

temperature  and  1  bar  (4.036  x  105  at  25°C) 
(mol  cm-3  bar-1) 

Ci  concentration  of  species  i  (mol  cm-3) 

c®  concentration  of  species  i  at  a  reference  condi¬ 

tion  (mol  cm-3) 

dCat  electrode  catalyst  loading  (mg  cm-2) 

Di  free  stream  diffusivity  of  species  i  (cm-2  s-1) 
D  f  effective  gas  diffusivity  of  species  i  (cm-2  s- 1 ) 

D'l  effective  liquid  phase  diffusivity  of  species  i 
(cm-2  s-1) 

£rev  theoretical  open  circuit  potential  at  the  given 
conditions  (V) 

E°  theoretical  open  circuit  potential  at  standard 
conditions  (V) 

F  Faraday’s  constant  (96485.309  C  mol-1) 

Hi  Henry’s  law  constant  for  species  i  (mol  cm-3 
bar-1) 

i  local  current  density  (A  cm-2) 

z'o  exchange  current  density  (A  cm-2) 

I  total  current  density  of  single  cell  (A  cm-2) 

Li  length  of  layer  i  (cm) 

mi  reaction  order  of  species  i 

Mi  symbol  for  the  chemical  formula  of  species  i 

n  number  of  electrons  transferred 

N i  molar  flux  of  species  i  (mol  cm-2  s-1) 

Pi  partial  pressure  of  species  i  (bar) 

p]  vapour  pressure  of  species  i  (bar) 

Q  ox  quantity  of  charge  to  reduce  surface  oxide 
(mC  cm-2  and  p,C  cm-2) 

R  ideal  gas  constant  (8.31451  JK-1  mol-1) 

R?  electrochemical  reaction  rate  per  unit  volume 
species  i  (mol  cm-3  s-1) 

R?  mass  transfer  rate  over  phase  boundary  species 
i  (mol  cm-3  s-1) 

Si  stoichiometric  coefficient  of  species  i 

t  time  (s) 

T  absolute  temperature  (K) 

ill  free  stream  mobility  of  species  i  (mol  cm2  J-1 

s-1) 

u'-  effective  mobility  of  species  i  (mol  cm2  J-1 
s-1) 

U  electrical  potential  of  an  electronically  conduc¬ 

tive  solid  phase  (V) 


f^cell 

cell  voltage  (V) 

V 

volume  average  velocity  (cms-1) 

Kat 

specific  catalyst  volume  (cm) 

Zi 

charge  number  of  species  i 

Greek 

letters 

aa 

anodic  transfer  coefficient 

a  c 

cathodic  transfer  coefficient 

8 

thickness  of  electrolyte  film  (cm) 

€ 

porosity 

T] 

local  overpotential  (V) 

Kc 

electronic  conductivity  of  cathode  (S  cm-1) 

Km 

ionic  conductivity  of  membrane  (S  cm-1) 

£ 

dimensionless  spatial  coordinate 

0 

solution  phase  potential  (V) 

Superscripts 

e 

electrochemical  reaction 

g 

gas  phase 

1 

liquid  phase 

P 

mass  transport  across  phase  boundary 

V 

vapour  pressure 

0 

reference  condition 

Subscripts 

(ads) 

adsorbed  specie 

(aq) 

solvated 

AN 

anode 

c 

cathode 

cat 

catalyst  (RUO2) 

CCL 

cathode  catalyst  layer 

DIFF 

cathode  diffusion  layer 

i 

species  i 

j 

species  j 

(1) 

liquid  phase 

m 

membrane 

MEM 

membrane 

SEP 

separator 

+ 

cations  (H+) 

anions  (Cl-) 

2.  Description  of  system 

The  H2-CI2-FC  single  cell  considered  in  this  study  con¬ 
sists  of  six  layers:  an  anode  gas  diffusion  layer,  an  anode 
catalyst  layer,  a  Nafion®  117  membrane,  a  separator  layer, 
a  cathode  catalyst  layer  and  a  cathode  gas  diffusion  layer.  A 
schematic  diagram  of  the  cell  is  shown  in  Fig.  1.  The  anode 
is  a  conventional  PEMFC-electrode,  based  either  on  a  carbon 
cloth  or  carbon  paper,  with  a  thin  layer  of  Pt/C  electrocatalyst 
bonded  to  a  Nafion®  117  membrane.  The  separator  layer  con¬ 
sists  of  a  porous  polyphenylene  sulphide  (PPS)  cloth  and  a 
circulating  hydrochloric  acid  electrolyte.  The  cathode  is  con- 
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Thomas- Alyea  [18].  It  is  also  assumed  that  in  the  gas  diffu¬ 
sion  electrode,  the  micropores  are  occupied  only  by  the  liquid 
phase  and  the  macropores  only  by  the  gas  phase. 

3.  Mathematical  modeling 

3.1.  Governing  equations 

The  equation  of  continuity  for  species  i  in  a  porous 
medium  can  be  written  in  the  general  form: 

-f-  =  —  VN;  +  Rf  +  (2) 


Fig.  1.  Schematic  diagram  of  the  H2-CI2-FC  single  cell. 

structed  of  a  highly  Teflon-impregnated  carbon  paper  acting 
as  the  gas  diffusion  layer  and  a  thin,  porous  catalytic  layer 
consisting  of  RuC>2  particles. 

The  most  distinct  feature  of  the  H2-CI2  fuel  cell  is  the 
cathode  catalyst  layer  and  the  separator  (electrolyte)  layer. 
To  simplify  the  model,  the  anode  gas  diffusion  and  catalyst 
layers  are  considered  to  play  a  negligible  role  in  the  overall 
performance  of  the  cell,  due  to  the  rapid  oxidation  kinetics 
and  high  diffusivity  of  hydrogen,  and  are  thus  omitted  from 
the  mathematical  model.  The  presence  of  the  liquid  elec¬ 
trolyte  and  its  intimate  contact  with  the  Nation  membrane 
will  probably  ensure  a  higher  degree  of  membrane  humid¬ 
ification  than  normally  found  in  PEM-fuel  cells.  Thus,  the 
Nation®  117  membrane  is  considered  to  be  fully  humidified 
and  acting  as  an  ideal  proton  conductor  with  a  constant  con¬ 
ductivity.  In  addition,  as  a  consequence  of  the  forced  convec¬ 
tion  of  the  liquid  electrolyte,  a  rapid  mass  transport  within  is 
to  be  expected.  It  is  therefore  assumed  that  the  mass  transfer 
of  the  electrolyte  and  its  dissolved  species  across  the  mem¬ 
brane  layer  is  negligible  and  that  the  membrane  thus  acts  as 
an  impermeable  barrier.  The  gaseous  chlorine  supplied  to  the 
cathode  gas  chamber  diffuses  across  the  cathode  gas  diffu¬ 
sion  layer  and  the  catalyst  layer  through  macropores  in  the 
electrode.  The  chlorine  then  dissolves  in  the  electrolyte  in 
the  catalyst  layer,  diffuses  through  the  electrolyte  film  cover¬ 
ing  the  catalyst  particles  and  reacts  electrochemically  on  the 
catalyst  surface  according  to  the  following  reaction: 

Cl2(aq)  +  2e-  -  2C1-  (1) 


where  e  is  porosity,  t  the  time,  q  and  N;  refer  to  the  con¬ 
centration  and  molar  flux  of  species  i,  respectively,  and  R f 
and  R ?  indicate  the  electrochemical  reaction  rate  and  mass 
transfer  rate  of  species  i  per  unit  volume  of  the  electrode. 

3.2.  The  flux  expression 

The  flux  expression,  N*,  depends  on  whether  the  species 
exists  in  the  gas  or  liquid  phase.  It  is  assumed  that  only  chlo¬ 
rine  and  hydrogen  chloride  exist  in  the  gas  phase  and  that 
these  species  have  an  ideal  gas  behaviour.  Thus,  Fick’s  law 
for  binary  diffusion  can  be  used 

Nf  =  -DfcgVPi  (3) 

where  D  f  is  the  effective  gas  diffusivity,  cg  the  molar  con¬ 
centration  of  an  ideal  gas  at  the  current  temperature  and  1  bar 
and  pl  is  the  partial  pressure  of  species  i.  For  the  liquid  phase, 
it  is  assumed  that  the  flux  equation  for  dilute  solutions  can 
be  employed: 

N)  =  -D'l  Vet  -  Ziu\Fci  S/0  +  c{s  (4) 

where  F  is  Faraday’s  constant,  0  the  solution  phase  potential, 
v  the  volume  average  velocity  and  Df,  Zi  and  u\  are  the  ef¬ 
fective  liquid  phase  diffusivity,  charge  number  and  effective 
mobility  of  species  /,  respectively. 

The  effective  diffusivity  and  mobility  is  corrected  for 
porosity  and  tortuosity  by  applying  the  Bruggemann  correc¬ 
tion  [18]: 

D\  =  15  (5) 

u'i  =  ute 1,5  (6) 


The  electrons  needed  by  the  above  reaction  move  from  the 
current  collector,  across  the  gas  diffusion  and  catalyst  layers 
where  they  are  consumed  by  the  chlorine  reduction  reaction. 
The  product,  chloride  ions,  either  diffuses  through  the  micro¬ 
pores  into  the  bulk  electrolyte  (separator  layer)  or  recombines 
with  protons  in  the  electrolyte  and  evaporates  as  gaseous  HC1 
and  diffuses  through  the  gas  diffusion  layer  and  into  the  gas 
channels.  Each  layer  of  the  cell  is  assumed  to  be  a  super¬ 
position  of  two  or  more  continua.  This  assumption  is  based 
on  the  porous  electrode  model  presented  by  Newman  and 


where  D\  and  ui  are  the  free  stream  diffusivity  and  mobility 
of  species  i. 


3.2.1.  Interfacial  mass  transfer 

With  the  assumption  of  equilibrium  at  the  gas-electrolyte 
interface,  the  mass  transport  rate,  R?,  for  species  i  across  a 
phase  boundary  can  be  approximated  to: 


RV  =  -a%D] 


1  (  Pi  Ci 


(7) 


8 
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where  ag  is  the  specific  area  of  the  gas-electrolyte  interface, 
Hf  the  Henry’s  law  constant  of  species  i  and  8  is  the  thickness 
of  the  electrolyte  film. 

3.2.2.  Electrochemical  reaction  rate 

The  electrochemical  reaction  rate  of  species  i  per  unit  vol¬ 
ume,  R f,  can  be  represented  by: 


where  a1  is  the  specific  area  of  the  catalyst-electrolyte  inter¬ 
face,  n  denotes  the  number  of  electrons  transferred  while  i  is 
the  local  current  density  and  Si  is  the  stoichiometric  coeffi¬ 
cient  of  species  i. 

The  stoichiometric  coefficient  is  given  by  expressing  an 
electrochemical  reaction  in  the  form: 


3.3.  Membrane  layer 

There  are  no  electrochemical  reactions  in  this  layer  and 
there  are  no  gas,  liquid  or  electronic  conductive  phases 
present.  It  is  assumed  that  the  potential  drop  in  this  layer 
can  be  described  by  Ohm’s  law. 

/ 

V0mem  =  —  (14) 

where  I  is  the  total  current  of  the  single  cell  and  /cm  is  the 
conductivity  of  the  membrane.  Since  there  are  no  electro¬ 
chemical  reactions  occurring  in  this  layer,  the  current  density 
is  constant.  By  differentiating  Eq.  (14),  the  potential  drop  in 
the  membrane  can  be  expressed  as: 

V20mem  =  o  (15) 


ne~  (9) 

i 

where  Mi  is  a  symbol  for  the  chemical  formula  of  species  i. 

The  local  current  density,  i,  can  be  described  by  the  Butler- 
Volmer  equation: 


i  =  i  o 


n 


c‘ , 

-0  )  exP 


n 


exp 


aa  nFrj 
RT 


-acnFrj 

RT 


(10) 


3.4.  Separator  (electrolyte)  layer 


There  are  no  electrochemical  reactions  in  this  layer  and 
there  are  no  gas  phase  or  electrically  conductive  solid  phase 
present.  Four  species  exist  in  this  layer:  liquid  water,  dis¬ 
solved  chlorine  and  two  ions —  and  Cl^.  By  employ¬ 

ing  the  continuity  equation  and  the  transport  equations  to  all 
species,  except  water,  one  obtains: 

^  —  V  (c+v)  =  0  (16) 


V(Zy+Vc+)  +  v(^^c+V0 


V(D'_  Vc_)  +  V 


/  z~D’_F 

V  RT 


C-V0 


V(c_v)  =  0  (17) 


where  /o  is  the  exchange  current  density,  cf  and  m,-  corre¬ 
spond  to  the  concentration  of  species  i  at  a  reference  con¬ 
dition  and  the  reaction  order  of  species  i  and  aa  and  ac  de¬ 
note  the  apparent  anodic  and  cathodic  transfer  coefficients, 
respectively. 

The  local  overpotential,  r\ ,  is  given  by: 

r)  =  U  —  @  —  £rev  (11) 

where  U  is  the  electrical  potential  of  an  electronically  conduc¬ 
tive  solid  phase  and  £rev  is  the  reversible  electrode  potential 
at  the  given  concentrations  given  by: 


V(Z^l2  VcCi2)  -  V(cCi2v)  =  0  (18) 

Since  there  are  no  reactions  or  mass  transfers  over  phase 
boundaries  in  this  layer,  there  is  no  change  in  the  volume  av¬ 
erage  velocity  (neglecting  the  influence  of  the  density  change 
due  to  HC1  dissolution): 

Vv  =  0  (19) 

Due  to  the  porous  PPS  cloth,  which  the  electrolyte  has  to 
flow  through,  it  is  assumed  that  the  electrolyte  flow  in  the 
y-direction  can  be  described  by  a  plug  flow  characteristic. 


£rev  =  E[J 


—  Ts, 
nF  4^  Vc° 


(12) 


where  E°  is  the  theoretical  open  circuit  potential  at  standard 
conditions. 


3.2.3.  Electroneutrality 

In  the  liquid  phase,  the  distribution  of  the  charged  species 
must  obey  the  electroneutrality  condition: 

T2  ziCi  =  o 


3.5.  Catalyst  layer 

Assuming  that  there  are  no  electrochemical  reactions  in 
the  gas  phase  of  the  catalyst  layer,  gaseous  chlorine  dissolves 
into  the  electrolyte  before  reacting.  The  product,  HC1,  either 
diffuses  into  the  bulk  electrolyte  or  evaporates  into  the  gas 
phase  and  diffuses  to  the  gas  channels.  Because  of  the  low 
vapour  pressure  of  water  at  25  °C,  it  is  assumed  that  there  is 
negligible  amounts  of  water  vapour  in  the  gas  phase;  thus, 
the  transport  equation  for  chlorine  gas  becomes: 


(13) 
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and  the  following  equation  can  be  obtained  from  the  conti¬ 
nuity  equation: 

V(£cVg  VpC|2)  -  alD'f  =  0  (21) 

In  the  liquid  phase,  the  mass  transport  of  the  dissolved  reac¬ 
tant  gases  can  be  written  by  using  the  Nernst-Planck  equation 
(Eq.  (4)).  Since  the  dissolved  chlorine  gas  is  an  electrically 
neutral  species,  the  second  term  on  the  right-hand  side  is  zero. 
The  dissolved  gas  reacts  electrochemically  at  the  catalyst- 
electrolyte  interface;  thus,  the  following  equation  is  obtained 
at  steady  state: 

V(^ci2  Vcci2)  -  V(cCi2v)  +  alD'f  ~  ^ 

_  ^C124;’c  =  Q 
ncF 

The  local  current  density  at  the  cathode,  ic,  is  described  by  ap¬ 
plying  the  Butler- Volmer  equation  to  the  chlorine  reduction 
reaction: 


The  transport  equation  for  gaseous  hydrochloric 
comes: 


(23) 

acid  be- 


NHCl  =  “^HClcg  vPhci 


(24) 


Thus,  with  the  interfacial  mass  transport  rate  of  HC1  being: 

*HC1  =  -«lAVg  (P"0scPH0)  (25) 

where  p^cl  is  the  concentration  dependent  vapour  pressure 
of  HC1,  the  following  equation  can  be  obtained  from  the  con¬ 
tinuity  equation: 

V^dCg  V^hci)  +  al D'£acg  (  P'HO~PHC'  )j  =  0  (26) 

Protons  do  not  react  electrochemically  at  the  cathode,  while 
the  electrochemical  reaction  rate  for  chloride  ions  follows 
Eq.  (8).  Thus,  using  the  continuity,  flux  and  Butler-Volmer 
equations,  we  get: 


V(Z)/_J_  Vc+)  +  z+F  V(w+c+V0)  -  V(c+v) 

-  «ogOHVg  (^C1~j?HC1)  =  0  (27) 


V  ( £>1  Vc_)  +  z-FV  (U-C-V&)  —  V  (c_v)  —  S—155l 
\  /  ncF 

~  «§  AVg  (^C1~/HC1)  =  0  (28) 

The  total  current  density,  /,  obtained  from  the  single  cell  must 
be  equal  to  the  integrated  value  of  the  local  current  density,  i, 
with  respect  to  the  thickness  (x)  of  the  cathode  catalyst  layer: 

fLcat 

/  =  -/  fli  dz  (29) 

Jo 

As  for  the  membrane  layer,  Ohm’s  law  can  describe  the  po¬ 
tential  drop  in  the  solid  catalyst  particles.  Differentiating  Eq. 
(29)  and  combining  it  with  Ohm’s  law  gives  the  following  ex¬ 
pression  for  the  potential  drop  in  the  cathode  catalyst  layer: 

V2t/  =  (30) 

where  /cc  is  the  effective  electrical  conductivity  of  the  cath¬ 
ode. 

3.6.  Gas  diffusion  layer 

There  are  no  electrochemical  reactions  in  this  layer; 
gaseous  chlorine  diffuses  from  the  gas  channels  towards  the 
catalyst  layer  while  the  product,  HC1,  diffuses  from  the  cata¬ 
lyst  layer  towards  the  gas  channels.  As  for  the  catalyst  layer, 
it  is  assumed  that  no  water  vapour  exists  in  the  gas  phase.  The 
transport  equation  for  chlorine  gas  and  gaseous  hydrochloric 
acid  becomes: 

Ni  =  —  D-cg  V Pi  (31) 

As  for  the  catalyst  layer,  the  potential  drop  can  be  described 
by  Ohm’s  law: 

S7U  =  —  (32) 

where  \c&  is  the  effective  conductivity  of  the  diffusion  layer. 
Since  no  electrochemical  reactions  occur  in  the  gas  diffusion 
layer,  I  is  constant.  By  differentiating  Eq.  (32),  the  potential 
drop  in  the  layer  can  be  expressed  as: 

V2l/diff  =  o  (33) 

3. 7.  Boundary  conditions 

3.7.1.  Anode/membrane  interface 

At  this  interface,  only  protons  are  transferred  and  the  flux 
of  protons  must  be  equal  to  the  total  current  density: 

N+|sep  =  —  (34) 

nF 

Since  the  anode  has  been  omitted  from  this  model,  the  elec¬ 
trical  potential  at  the  anode/membrane  interface  can  be  arbi- 
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trarily  fixed.  Hence,  the  membrane  potential  at  this  boundary 
is  set  to  zero: 

^mem  I  AN  =  0  (35) 

3.7.2.  Membrane/ separator  interface 

The  membrane  acts  as  an  impermeable  barrier  for  the  elec¬ 
trolyte  and  its  dissolved  species.  Thus,  the  fluxes  of  anions 
and  chlorine  are  zero,  while  the  flux  of  protons  is  continuous 
at  this  interface 

Nci2(aq)lsEP  =  o  (36) 

N-|Sep  =  0  (37) 

N+|mem  =  N+|sep  (38) 

The  solution  potential  is  equal  to  the  membrane  potential, 
^mem  I  MEM  =  ^sepIsEP  (39) 


3.7.3.  Separator/cathode  catalyst  interface 

It  is  assumed  that  the  PPS-separator  prevents  gaseous  re¬ 
actants  from  flowing  towards  the  opposite  electrode.  Thus, 
the  fluxes  of  gaseous  chlorine  and  hydrochloric  acid  are  zero 


at  this  interface: 

N;(g)lccL  =  0  =>  V/?|  Iccl  =  0  (40) 

Since  the  separator  is  an  electronic  insulator,  the  electronic 
current  density  in  this  layer  is  zero: 

V  U  |  ccl  =  o  (41) 

The  fluxes  of  dissolved  chlorine  and  ions  are  continuous  at 
this  interface: 

N+|ccl  =  N+|sep  (42) 

N-|ccl  =  N_|sep  (43) 

Nci2(aq)  I  ccl  =  Nci2(aq)  I  sep  (44) 


3.7.4.  Cathode  catalyst/gas  diffusion  layer  interface 

Since  the  gas  diffusion  layer  prevents  leakage  of  elec¬ 
trolyte,  it  is  assumed  that  the  fluxes  of  the  electrolyte  and 
its  dissolved  species  are  zero  at  this  boundary.  The  fluxes  of 
chlorine  gas  and  gaseous  hydrochloric  acid  are  continuous  at 


this  boundary. 

N/(g)  Idiff  =  N;(g)  Iccl  (45) 

N+|diff  =  0  (46) 

N_  Idiff  =  0  (47) 

Na2(aq)  Idiff  =  0  (48) 


3.7.5.  Gas  diffusion  layer/gas  channel  interface 

The  electrical  potential  can  be  arbitrarily  fixed  either  at 
the  anode  or  cathode  current  collector.  Setting  the  value  of 
the  anode  to  zero,  Uc  at  this  interface  will  be: 


Uq  —  f^cell 


(49) 


where  £/ce n  is  the  total  cell  voltage. 

Assuming  a  uniform  distribution  of  gases  in  the  gas  chan¬ 
nels,  the  partial  pressure  of  chlorine  and  hydrochloric  acid 
can  be  set  to: 


PC12  =  Pc  12  (5°) 

Phci  =  0  (51) 


3.7.6.  Electrolyte  inlet 

At  y  =  0,  the  electrolyte  concentration  and  volume  aver¬ 
age  velocity  are  known.  It  is  also  assumed  that  the  chlorine 
concentration  in  the  supplied  electrolyte  is  equal  to  1%  of 
the  saturation  amount  since  the  electrolyte  is  recycled  and 
removal  of  all  chlorine  gas  from  this  electrolyte  is  difficult. 

Cely=o  =  c°e  (52) 

v*l;y=o  =  v$  (53) 

Cc\2\y=0  =  0.01  •  Cci2,sat  (54) 


3. 7. 7.  Electrolyte  outlet 

At  y  =  1 ,  the  change  in  electrolyte  and  chlorine  concen¬ 
tration  is  zero.  It  is  assumed  that  the  mass  transport  due  to 
diffusion  and  migration  is  negligible  compared  to  the  con¬ 
vective  flux. 

V  (ce  |  y=  0  =  0  (55) 

V(CCl2(aq)l;y=l)  =  0  (56) 

W(\y\y=i)  =  0  (57) 


4.  Model  parameters  and  correlation 

Using  realistic  values  for  the  modeling  parameters  is  es¬ 
sential  for  the  viability  of  the  mathematical  model.  A  reason¬ 
able  choice  of  the  parameter  values  is  thus  as  important  as  the 
mathematical  modeling  itself.  Hence,  an  effort  has  been  made 
to  supply  as  realistic  values  as  possible  for  these  parameters. 

4.1.  Thermodynamic  properties 

4.1.1.  Solubility  of  chlorine 

It  is  assumed  that  the  solubility  of  chlorine  gas  obeys 
Henry’s  law: 

et  =  HiPi  (58) 

The  solubility  of  chlorine  changes  linearly  with  the  hy¬ 
drochloric  acid  concentration  in  the  range  applicable  for  this 
model.  Using  the  data  from  [19],  the  following  correlation 
was  obtained  for  the  solubility  of  chlorine  in  hydrochloric 
acid  at  25  °C: 

7/q2(^h+)  =  0-0027  •  cH+  +  6.0  x  10  ^ 


(59) 
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4.1.2.  Vapour  pressure  of  HCl 

The  vapour  pressure  of  HCl  over  hydrochloric  acid  in¬ 
creases  exponentially  with  increasing  hydrochloric  acid  con¬ 
centration.  At  25  °C,  the  following  correlation  was  obtained 
using  the  data  from  Perry ’s  Chemical  Engineers  ’  Handbook 
[20]: 

Phci(ch+)  =  2-4  x  10“6  exp(961.46  •  cH+)  (60) 

4.1.3.  Transport  properties 

The  mobility  of  the  electrolyte  species  changes  with  elec¬ 
trolyte  concentration.  Using  the  data  in  [21],  the  following 
correlations  were  obtained: 

ucr  =  7  X  1CT9  exp(— 171.51  •  ce)  (61) 

uu+  =  4  x  10“8  exp(— 171.51  •  ce)  (62) 

The  diffusivities  of  the  charged  species  were  found  using  the 
Nemst-Einstein  equation: 

Di  =  RTm  (63) 


Table  1 


Electrochemical  kinetic  parameters 


Parameter 

Value 

Number  of  electrons  transferred,  n 

2 

Stoichiometric  coefficient  of  chlorine,  sq\2 

-1 

Stoichiometric  coefficient  of  chloride,  scl- 

2 

Reaction  order  of  chlorine,  mci2 

1 

Reaction  order  of  chloride,  mcl- 

2 

Anodic  transfer  coefficient,  aa 

0.75 

Cathodic  transfer  coefficient,  ac 

0.25 

Exchange  current  density,  io  (A  cm-2) 

1.0  x  IO-5 

4.1.5.  Structural  parameters 

The  thickness  of  the  cathode  catalyst  layer  was  found  to 
be  10  p,m  by  analysing  a  cross-section  of  an  electrode  in  a 
scanning  electron  microscope.  The  gas  diffusion  electrode, 
Toray  carbon  paper  had  a  thickness  of  0.4  mm.  The  specific 
surface  area  of  the  Ru02  catalyst  was  45-65  m2  g- 1  (Alpha 
Aesar)  and  the  catalyst  loading  was  2  mg  cm-2.  By  assuming 
perfect  wetting  of  the  whole  catalyst  surface,  the  specific 
catalyst-electrode  interface  area,  a1,  can  be  calculated  from: 


The  diffusion  coefficient  for  chlorine  is  assumed  not  to  be 
influenced  by  the  acid  concentration.  The  diffusion  coeffi¬ 
cient  in  0.1N  HCl  at  25  °C  has  been  found  to  be  1.38  x 
10-5  cm2  s-1  [22].  The  electrical  conductivities  of  the  cata¬ 
lyst  and  gas  diffusion  layers  are  reported  to  be  7.14  S  cm-1 
[23]  and  12.5  S  cm-1,  respectively.  Earlier  work  from  our 
group  [24]  investigated  the  conductivity  of  Nation®  117  and 
found  a  value  of  0.1  S  cm-1  at  25  °C  in  a  0.3  mol  dm-3  HCl 
solution. 

4.1.4.  Kinetic  parameters 

The  exchange  current  density,  io,  for  chlorine  reduction  on 
Ru02  has  been  measured  to  be  ~1.0x  10-5  in  1  mol  dm-3 
HCl  saturated  with  Cl2  using  the  rotating  disc  electrode  tech¬ 
nique  [25].  The  anodic  and  cathodic  Tafel  slopes  for  the  chlo¬ 
rine  evolution/reduction  reaction  on  Ru02  have  been  reported 
by  several  authors  to  be  40  and  120  mV  dec-1,  respectively 
[26,27].  Thus,  the  anodic  and  cathodic  transfer  coefficients, 
ckfa  and  ac ,  are  0.75  and  0.25,  respectively.  Since  the  electro¬ 
chemical  reactions  within  the  cathode  catalyst  layer  are  con¬ 
sidered  to  be  elementary  reactions,  the  same  absolute  values 
as  the  stoichiometric  coefficients  are  used  as  reaction  orders 
of  the  reactants  and  products.  The  electrochemical  kinetic 
parameters  are  listed  in  Table  1 . 


where  JCat  is  the  catalyst  loading  (g  cm-2),  Vcat  the  specific 
catalyst  volume  (cm)  and  Acat  is  the  catalyst  surface  area 
(cm2  g-1). 

Hence,  the  value  of  a1  in  the  cathode  catalyst  layer  has 
a  value  of  9  x  105  cm2  cm-3.  Several  researchers  have  re¬ 
ported  specific  gas-electrolyte  interface  areas  of  the  order  of 
103  cm2  cm-3  and  an  electrolyte  film  thickness  of  ~0.5  p,m 
[21,28].  Thus,  the  values  for  ag  and  8  have  been  set  to 
1.0  x  103  and  0.5  p,m,  respectively.  The  structural  param¬ 
eters  are  listed  in  Table  2. 

4.1.6.  Operating  conditions 

The  operating  conditions  for  the  laboratory  test  cell  are 
ambient  temperature  and  pressure.  The  base  case  operating 
conditions  for  the  model  are  listed  in  Table  3. 

4.2.  Method  of  solution 

The  model  equations  are  highly  coupled  and  non-linear 
and  a  numerical  solution  is  required.  The  finite  element 
method  program  Femlab®  3.0a  with  the  chemical  engineer- 


Table  2 

Base  case  structural  parameters 


Parameter  (unit) 

Membrane  layer 

Separator  layer 

Catalyst  layer 

Gas  diffusion  layer 

Layer  thickness,  Li  (cm) 

0.0175 

0.08 

1.0  X  10“3 

0.04 

Gas  phase  porosity,  eg 

0.1 

0.7 

Liquid  phase  porosity,  d 

0.5 

0.3 

Solid  phase  porosity,  es 

0.6 

0.3 

Specific  area  of  gas-electrolyte  interface,  ag  (cm-1) 

1.0  x  103 

Specific  area  of  catalyst-electrolyte  interface,  a 1  (cm-1) 

1.0  x  105 

Thickness  of  electrolyte  film,  8  (cm) 

5.0  x  10“4 

Electrical  conductivity  (S  cm-1) 

0.10 

7.14 

12.5 
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Table  3 

Base  case  operating  conditions 

Parameter  (unit) 

Value 

Electrolyte  inlet  concentration,  Ce  (mol  cm-3) 

0.003 

Operating  temperature,  T  (°C) 

25 

Operating  pressure,  p°  (bar) 

1 

Electrolyte  flow  rate,  v  (cms-1) 

0.02 

in g  toolbox  from  Comsol  AB  was  employed  to  solve  the 
equations.  A  multiphysics  model  incorporating  three  appli¬ 
cation  modes  and  six  dependent  variables  was  created  and 
a  weak,  non-linear  parametric  solver  was  used  to  solve  the 
constructed  model. 


5.  Experimental  validation 

The  experimental  single  cell  is  a  further  development  of 
a  cell  examined  in  [29].  Fig.  2  shows  a  schematic  illustra¬ 
tion  of  the  fuel  cell  construction.  The  housing  consisted  of 
two  graphite  current  collectors  with  double  serpentine  flow 
fields,  a  carbon  paper  anode  attached  to  a  Nation  117  mem¬ 
brane,  a  polyphenylene  sulphide  (PPS)  cloth  separator,  a  car¬ 
bon  paper  cathode  and  a  silicon  gasket.  The  graphite  current 
collectors  included  inlets  and  outlets  for  the  circulating  liquid 
HC1  electrolyte.  The  PPS  cloth  acts  as  a  spacer  between  the 
electrodes  making  it  possible  to  apply  a  mechanical  pressure 
to  the  electrodes  and  simultaneously  supplying  space  for  the 
circulating  electrolyte. 

5.7.  Cathode  preparation 

The  cathodes  were  constructed  of  ELAT  carbon  paper 
(E-TEK)  with  various  amounts  of  Teflon  impregnation  (0- 


Nafion 


Fig.  2.  Schematic  illustration  of  the  hybrid  fuel  cell  construction. 


60  wt%).  Catalyst  inks  were  prepared  using  commercial  un¬ 
supported  Ru02  catalyst  (Alpha  Aesar).  The  catalyst  was 
mixed  with  various  amounts  of  Teflon  suspension  and  iso¬ 
propanol  to  produce  a  range  of  catalysts  with  different  degrees 
of  wet-proofing.  The  mixture  was  stirred  and  ultrasonicated 
for  several  hours  before  application  to  the  cathode  carbon 
paper  electrode  using  an  air-brush  (Badger),  with  the  carbon 
electrode  placed  on  a  thermal-controlled  heating  plate.  After 
spraying,  the  electrode  was  dried  in  a  forced  convection  oven 
at  150  °C  for  1  h. 

5.2.  Membrane  and  ME  A  preparation 

The  Nation  117  membranes  were  treated  by  immersion  in 
85  °C  18.2  distilled  water  for  15  min,  thereafter  in  5% 
H2O2  (Merck  p. a.)  for  30  min  and  washed  with  distilled  wa¬ 
ter  and  ion  exchanged  twice  in  0.05  M  H2SO4  (Merck  p.a.), 
each  time  for  30  min.  After  ion  exchange,  the  membranes 
were  washed  thoroughly  by  immersion  in  distilled  water  four 
times,  each  for  about  15  min.  The  membranes  were  then  cut 
in  rectangular  pieces  and  stored  in  purified  water.  Commer¬ 
cially  available  ELAT  gas  diffusion  electrodes  (E-TEK)  with 
a  20  wt%  Pt/C  catalyst  loading  of  1  mg  Pt  cm-2  were  used  as 
anodes.  A  rectangular  piece  of  Nation  membrane  was  gen¬ 
tly  dried  with  a  lint  free  paper,  sandwiched  with  a  6.5  cm2 
rectangular  anode  piece  between  two  thin  Teflon  sheets  and 
hot-pressed  for  3  min  at  130  °C  and  10  MPa.  The  MEA  was 
then  cooled  down  and  placed  in  purified  water  until  mounted 
in  the  cell. 

5.3.  Fuel  cell  testing 

The  cell  was  assembled  using  a  torque  of  0.15  Nm  on 
the  bolts.  All  experiments  reported  here  were  conducted  at 
room  temperature  and  atmospheric  pressure.  The  gases,  hy¬ 
drogen  (5.0  AG  A)  and  chlorine  (2.8  Gerling  Holtz),  were  fed 
to  the  cell  without  any  kind  of  pre-heating  or  humidification 
and  exited  the  cell  through  a  water  column  to  create  a  small 
overpressure  inside  the  cell  (10  cm  H20).  The  fuel  cell  was 
controlled  using  a  Solartron  SI  1287  electrochemical  inter¬ 
face.  Quasi- steady- state  polarization  curves  were  recorded  at 
an  operating  temperature  of  25  °C  by  sweeping  the  potential 
at  lmV  s-1  from  open  circuit  to  approximately  0.2  V  and 
then  reversed.  The  HC1  liquid  electrolyte  circulated  from  a 
thermally  controlled  external  reservoir  through  the  cell  at  a 
rate  of  0.2  ml  s-1  via  a  peristaltic  pump.  The  HC1  electrolyte 
concentration  was  3  mol  dm-3. 


6.  Results  and  discussion 

The  model  developed  is  a  two-dimensional  model  and  the 
fuel  cell  performance  and  the  distribution  of  the  physical  pa¬ 
rameters  will  vary  along  the  cell  length  (y-axis).  It  was  found 
that,  for  the  base  case  electrolyte  velocity,  this  variation  was 
negligible  from  a  distance  of  0.05  cm  from  the  electrolyte 
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Fig.  3.  Polarization  curve  for  base  case.  The  open  squares  are  experimental 
results  from  a  6  cm2  laboratory  cell.  The  cathodic  overpotential  (rjc)  and  iR 
drop  of  the  cell  are  also  shown. 

inlet  for  all  current  densities.  The  results  presented  in  this 
section  are  taken  from  a  point  0.5  cm  from  the  electrolyte 
inlet  and  are  considered  to  give  a  good  representation  of  the 
distribution  of  the  physical  parameters  inside  the  whole  cell. 

6.1.  Comparison  of  model  with  experimental  results 

A  comparison  of  the  model  prediction  for  the  base  case 
with  experimental  data  obtained  with  the  laboratory  test  cell 
described  above  is  presented  in  Fig.  3.  The  open  squares  rep¬ 
resent  the  mean  values  from  five  experimental  runs  with  er¬ 
ror  bars  giving  the  corresponding  deviation  from  the  mean 
values,  while  the  full  line  is  the  model  prediction.  The  rela¬ 
tively  large  uncertainty  in  the  experimental  results  is  due  to  a 
very  unstable  operation  and  low  experimental  reproducibil¬ 
ity  of  the  laboratory  cell.  The  experimental  aspects  of  this 
cell  will  be  discussed  in  more  detail  in  an  upcoming  paper. 
Furthermore,  no  experimental  results  were  recorded  above 
0.35  A  cm-2  due  the  potentiostat  limiting  current  range  of 

2  A. 

6.2.  Polarization  curve  for  base  case 

The  polarization  curve  for  the  base  case  presented  in  Fig. 

3  displays  a  typical  shape  for  low  temperature  fuel  cells.  It 
includes  two  of  the  three  distinguishable  zones:  activation 
controlled,  ohmic  controlled  and  mass  transport  controlled 
polarization.  The  activation  controlled  zone  is  much  less  pro¬ 
nounced  than  in  similar  polarization  curves  for  H2-O2  fuel 
cells.  This  is  most  probably  because  the  kinetics  of  chlo¬ 
rine  reduction  is  faster  than  the  kinetics  for  oxygen  reduc¬ 
tion;  the  exchange  current  density  for  chlorine  reduction  on 
ruthenium  oxide  is  ~1.0x  10-5Acm-2,  compared  to  the 
exchange  current  density  of  oxygen  reduction  on  platinum, 
which  has  a  value  of  ~  1 .0  x  10-9  A  cm-2  [30] .  The  slope  of 
the  linear  region  corresponds  to  a  cell  resistance  of  approxi¬ 
mately  0.56  Q  cm2.  The  fuel  cell  has  an  open  circuit  voltage 


of  1.26  V  and  at  a  cell  voltage  of  0.2  V,  no  limiting  current 
can  be  observed.  Although  a  slight  increase  of  the  slope  of  the 
polarization  curve  can  be  seen.  The  maximum  power  density 
of  the  cell  is  about  0.58  W  cm-2  at  a  cell  voltage  of  0.55  V 
(determined  from  the  mathematical  model).  The  iR  drop  in 
the  membrane  and  electrolyte  and  the  cathodic  overpoten¬ 
tial  predicted  by  the  mathematical  model  are  also  shown  in 
Fig.  3.  The  overpotential  dominates  the  losses  at  low  current 
densities,  while  at  current  densities  above  0.6  A  cm-2,  the 
iR  drop  is  responsible  for  more  than  50%  of  the  polarization 
losses.  At  very  high  current  densities  (>1 .6  A  cm-2),  the  cath¬ 
ode  overvoltage  starts  to  increase  more  rapidly,  indicating  an 
impending  mass  transport  controlled  limiting  current. 

6.3.  Distribution  of  physical  variables 

The  distribution  of  the  partial  pressure  of  gaseous  chlo¬ 
rine  as  a  function  of  cell  voltage  and  dimensionless  spatial 
coordinate,  §,  is  given  in  Fig.  4.  It  can  be  seen  that  with 
decreasing  cell  voltage,  the  distribution  of  the  partial  pres¬ 
sure  of  chlorine  becomes  less  uniform.  At  a  cell  voltage  of 
0.2  V,  the  partial  pressure  of  chlorine  inside  the  catalyst  layer 
is  approximately  0.85  atm  with  a  small  gradient  towards  the 
diffusion  layer.  This  indicates  that  even  though  there  exists  a 
certain  mass  transfer  resistance  for  the  chlorine  gas,  it  is  not 
the  factor  limiting  the  total  cell  performance.  The  distribu¬ 
tion  of  dissolved  chlorine  is  shown  in  Fig.  5.  At  cell  voltages 
between  1.2  and  0.9  V,  there  is  only  a  very  slight  decrease 
of  the  chlorine  concentration.  With  a  further  lowering  of  the 
cell  voltage,  an  interesting  phenomenon  occurs:  the  concen¬ 
tration  of  dissolved  chlorine  increases  and  reaches  a  peak 
value  at  Uce n  ~  0.75  V.  Below  0.75  V,  the  chlorine  concen¬ 
tration  rapidly  decreases.  This  unexpected  behaviour  is  due 
to  the  increasing  solubility  of  chlorine  with  increasing  elec¬ 
trolyte  concentration  (see  Fig.  6  and  Eq.  (59)).  At  0.75  V,  the 
electrolyte  concentration  reaches  a  plateau  value  and  low- 

Cathode  Diffusion 


Fig.  4.  Chlorine  partial  pressure  distribution  in  the  cathode  catalyst  and 
diffusion  layer  as  a  function  of  cell  voltage. 
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Cathode  Catalyst 


Fig.  5.  Dissolved  chlorine  concentration  distribution  in  electrolyte  and  cath¬ 
ode  catalyst  layers  as  a  function  of  cell  voltage. 

ering  the  cell  voltage  further  increases  the  electrochemical 
reduction  rate  of  the  chlorine  reduction  reaction  and  the  con¬ 
centration  of  dissolved  chlorine  drops.  The  low  concentration 
of  dissolved  chlorine  at  low  cell  voltages  indicates  that  this 
parameter  is  the  cause  of  the  tendency  of  a  limiting  current 
behaviour  observed  in  Fig.  3.  Fig.  6  displays  the  concentra¬ 
tion  distribution  of  HC1  in  the  electrolyte  as  a  function  of 
cell  voltage.  It  can  be  seen  that  there  is  a  drastic  increase  in 
acid  concentration  with  decreasing  cell  voltage  down  to  0.6  V 
where  it  reaches  a  plateau  value.  This  increase  is  caused  by 
the  production  of  chloride  ions  from  the  reduction  of  chlo¬ 
rine  and  the  diffusional  and  migrational  fluxes  of  chloride 
ions  are  too  low  to  maintain  a  concentration  close  to  the  in¬ 
let  concentration.  This  high  electrolyte  concentration  inside 
the  catalyst  layer  has  several  effects:  it  causes  a  drop  in  the 
reversible  voltage  of  the  cathodic  reaction  and  exposes  the 


Cathode 


Fig.  6.  Electrolyte  concentration  distribution  in  the  electrolyte  and  cathode 
layer  as  a  function  of  cell  voltage. 


Cathode  Catalyst 


Fig.  7.  HC1  partial  pressure  distribution  in  the  cathode  catalyst  and  diffusion 
layer  as  a  function  of  cell  voltage. 

cathode  catalyst  to  a  more  corrosive  environment.  Another 
effect,  not  taken  into  consideration  in  this  work,  is  the  pos¬ 
sibility  of  increased  adsorption  of  Cl-  with  increasing  HC1 
concentration  which  may  reduce  the  chlorine  reduction  rate. 
The  electrolyte  has  a  maximum  conductivity  at  concentra¬ 
tions  close  to  5  mol  dm-3;  thus,  the  concentration  build-up 
inside  the  catalyst  layer  has  a  negative  effect  on  the  conduc¬ 
tivity  in  this  layer,  causing  a  higher  ohmic  drop. 

On  the  other  hand,  a  positive  effect  of  the  concentration 
increase  is  a  higher  solubility  of  chlorine  in  the  catalyst  layer. 
Fig.  7  shows  the  distribution  of  the  partial  pressure  of  gaseous 
HC1  in  the  catalyst  and  diffusion  layers  as  a  function  of  cell 
voltage.  The  evaporation  of  HC1  has  a  strong  correlation  with 
the  concentration  of  hydrochloric  acid  in  the  catalyst  layer 
shown  in  Fig.  6.  This  phenomenon  manifests  itself  only  at 
low  cell  voltages,  where  the  partial  pressure  of  hydrochlo¬ 
ric  acid  reaches  values  of  approximately  0.15  atm.  It  does 
not  seem  that  the  mass  transport  resistance  of  gaseous  hy¬ 
drochloric  acid  out  of  the  catalyst  layer  has  any  major  effect 
on  the  cell  performance.  Figs.  8  and  9  show  the  profiles  of  the 
local  apparent  current  density  and  local  overvoltage,  respec¬ 
tively.  The  distribution  of  these  parameters  is  very  uniform 
through  the  catalytic  layer  at  all  cell  voltages,  it  is  only  at  low 
cell  voltages  that  a  small  curvature  close  to  the  separator  layer 
can  be  seen  on  the  local  current  density  distribution.  The  lo¬ 
cal  overvoltage  of  the  cathode  reaches  a  value  of  —0.36  V  at 
a  cell  voltage  of  0.2  V  and  does  not  have  the  typical  Tafelian 
behaviour.  This  is  probably  due  to  the  increasing  electrolyte 
concentration  and  the  depletion  of  dissolved  chlorine  in  the 
catalyst  layer  at  low  cell  voltages,  which  effectively  increases 
the  overvoltage  value  additionally.  Fig.  10  shows  the  profile 
of  the  membrane  and  solution  potentials  as  a  function  of  cell 
voltage.  The  variation  across  the  membrane  and  separator 
layers  corresponds  to  the  ohmic  drop  between  the  electrodes. 
The  potential  drop  in  the  cathode  catalyst  layer  is  negligi¬ 
ble  compared  to  the  quite  large  drops  in  the  separator  and 
membrane  layers.  The  figure  also  shows  that  although  the 
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Fig.  8.  Profile  of  the  local  current  density  in  the  cathode  layer  as  a  function 
of  cell  voltage. 


separator  layer  is  almost  five  times  thicker  than  the  mem¬ 
brane  layer,  the  potential  drop  across  the  two  layers  is  nearly 
equal. 

6.4.  Influence  of  operating  conditions 

6.4.1.  Initial  electrolyte  concentration 

Since  most  thermodynamic,  transport  and  kinetic  param¬ 
eters  are  influenced  by  the  electrolyte  concentration,  a  varia¬ 
tion  of  the  initial  electrolyte  concentration  can  have  a  strong 
effect  on  the  cell  performance.  For  instance,  increasing  the 
electrolyte  concentration  will  increase  the  solubility  of  chlo¬ 
rine,  but  above  5  mol  dm-3 ,  the  conductivity  will  decrease.  In 
addition,  an  increase  of  the  chloride  concentration  will  affect 
the  cathode  kinetics  by  increasing  the  influence  of  the  an¬ 
odic  back  reaction  and  reducing  the  theoretical  open  circuit 
voltage,  if  rev- 


Fig.  9.  Profile  of  the  local  overvoltage  in  the  cathode  layer  as  a  function  of 
cell  voltage. 
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Fig.  10.  Profile  of  the  membrane  and  solution  potentials  as  a  function  of  cell 
voltage. 

Polarization  curves  for  a  single  cell  with  different  elec¬ 
trolyte  inlet  concentrations  are  given  in  Fig.  11.  A  consider¬ 
able  increase  in  the  cell  performance  can  be  seen  when  in¬ 
creasing  the  electrolyte  concentration  from  1  to  3  mol  dm-3 
and  as  expected,  the  open  circuit  voltage  drops  somewhat. 
A  further  increase  from  3  to  5  mol  dm-3  has  an  negligible 
effect.  This  behaviour  was  observed  previously  on  a  fuel  cell 
similar  to  the  experimental  cell  presented  in  this  work  [29]. 
In  Fig.  12,  the  effect  of  initial  electrolyte  concentration  on 
current  density  as  a  function  of  cell  voltage  is  presented. 
At  high  cell  voltage,  the  variation  of  electrolyte  inlet  con¬ 
centration  has  a  very  low  effect  on  the  total  current  density 
of  the  cell.  At  lower  cell  voltages,  the  current  density  in¬ 
creases  significantly  with  increasing  concentration  from  1  to 
3  mol  dm-3  and  moves  through  a  broad  maximum  at  con¬ 
centrations  around  6  mol  dm-3.  The  dashed  line  in  Fig.  12 
shows  the  conductivity  of  hydrochloric  acid  as  a  function  of 
concentration.  It  is  clear  that  although  the  current  density  con¬ 
centration  dependence  has  the  same  trend  as  the  electrolyte 
conductivity,  the  effect  is  not  as  pronounced.  This  indicates 
that  other  mechanisms,  such  as  electrode  kinetics  and  mass 
transport,  influence  the  electrolyte  concentration  dependence 
of  the  current  density. 

6.4.2.  Operating  pressure 

Polarization  curves  for  the  single  cell  at  different  operat¬ 
ing  pressures  are  presented  in  Fig.  13.  An  increase  in  the  cell 
current  density  at  all  voltages  can  be  observed,  while  the  po¬ 
larization  curves  maintain  essentially  the  same  slope  in  the 
ohmic  polarization  region.  The  small  tendency  towards  lim¬ 
iting  current  behaviour  observed  at  an  operating  pressure  of 
1  bar  is  completely  eliminated  at  higher  pressures.  The  rea¬ 
son  for  the  increased  cell  performance  is  a  higher  interfacial 
mass  transport  of  chlorine  from  the  gas  phase  to  the  liquid 
phase  and  a  higher  solubility  of  chlorine  in  the  electrolyte. 
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Fig.  11.  Single  cell  polarization  curves  at  different  electrolyte  inlet  concentrations. 
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Fig.  12.  Cell  current  density  as  a  function  of  electrolyte  inlet  concentration  and  cell  voltage.  Dashed  line  represents  the  conductivity  of  HC1  as  a  function  of 
electrolyte  concentration. 


Fig.  13.  Single  cell  polarization  curves  at  different  operating  pressures. 
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There  is  no  indication  that  the  conductivity  of  the  electrolyte 
varies  with  the  operating  pressure. 

7.  Conclusions 

A  mathematical  model  for  a  H2-CI2  single  fuel  cell 
has  been  developed.  The  model  describes  mass  transport 
and  electrochemical  reactions  occurring  in  the  cathode  gas 
diffusion  layer,  the  cathode  catalyst  layer  and  the  separa¬ 
tor/membrane  layers.  From  the  profiles  of  dissolved  chlorine 
as  a  function  of  cell  voltage,  it  is  determined  that  the  inter¬ 
facial  mass  transport  of  chlorine  between  the  gas  and  liq¬ 
uid  phases  in  the  cathode  probably  is  the  rate-determining 
step  at  high  current  densities.  It  is  found  that  concentration 
of  the  HC1  electrolyte  inside  the  cathode  catalyst  layer  ap¬ 
proaches  values  close  to  that  of  concentrated  hydrochloric 
acid  (12  mol  dm-3),  especially  at  high  current  densities. 

Influences  of  the  operating  condition  of  the  fuel  cell  are 
investigated.  It  is  shown  that  the  cell  performance  increases 
substantially  by  increasing  the  inlet  concentration  of  the  hy¬ 
drochloric  acid  electrolyte  from  1  to  3  mol  dm-3.  Increas¬ 
ing  the  electrolyte  concentration  further  only  leads  to  small 
changes  in  the  overall  cell  performance.  By  increasing  the 
operating  pressure,  a  steady  increase  of  the  cell  performance 
due  to  faster  electrode  kinetics  can  be  observed. 
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